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Abstract 

A general theory of preparational uncertainty relations for a quanfum parficle in 
one spatial dimension is developed. We derive conditions which determine whether 
a given smooth function of fhe parficle's variances and ifs covariance is bounded 
from below. Whenever a global minimum exisfs, an uncerfainfy relafion has been 
obfained. The squeezed number sfafes of a harmonic oscillator are found to be univer¬ 
sal: no ofher pure or mixed sfafes will safurafe any such relafion. Geomefrically, we 
idenfify a convex uncertainty region in fhe space of second momenfs which is bounded 
by fhe inequalify derived by Roberfson and Schrodinger. Our approach provides a 
unified perspecfive on exisfing uncerfainfy relations for a single confinuous variable, 
and if leads to new inequalifies for second momenfs which can be checked experi¬ 
mentally. 


1 Introduction and main result 

Inspired by Heisenberg's analysis flj of Compton scattering, Kennard [2J proved the 
preparational uncertainty relation 

( 1 ) 

for the standard deviations Ap and Aq of momentum and position of a quantum particle 
with a single spatial degree of freedom. Experimentally, they are determined by mea¬ 
surements performed on an ensemble of systems prepared in a specific state \xp). The 
only states which saturate the bound Q are squeezed states with a real squeezing parame¬ 
ter ||3l|4l 13 (we follow the review lb] regarding the naming of squeezed states). Squeezed 
states are conceptually important since they achieve the best possible localization of a 
quantum particle in phase-space, and they are easily visualized by "ellipses of uncer¬ 
tainty". Each squeezed state may be displaced rigidly in phase space without affecting 
the value of the variances, resulting in a three-parameter family of states saturating the 
lower bound Q. 

Not many other uncertainty relations are known. The sum of the position and mo¬ 
mentum variances is bounded IZlISl according to the relation 
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A^p + A^qyfi, ( 2 ) 

which holds in a system of units where the physical dimensions of both position and 
momentum equal Vh. Only the ground state of a harmonic oscillator with unit mass and 
frequency saturates this inequality (ignoring rigid displacements in phase-space). The 
Robertson-Schrodinger (RS) inequality ll9l[T0l, 

A^pA^q-Cl^> (3) 

sharpens Heisenberg's inequality ([T} by including the covariance Cpq defined in Eq. Q. 
Eq.§ is saturated by the two-parameter family of squeezed states with a complex squeez¬ 
ing parameter |j6|, again ignoring phase-space displacements. The additional free pa¬ 
rameter describes the phase-space orientation of the uncertainty ellipse which, in the 
previous case, was aligned with the position and momentum axes. 

By introducing the observable f = —p — q, which satisfies the commutation relations 
[q, f] = [f,p] =h/i, one obtains a bound on the product of the variances of three pairwise 
canonical observables, 

A^p A^q A^r > ^ (x) ^ 

This triple product uncertainty relation has been found only recently Ifl^ . Ignoring phase- 
space translations, only one state exists which achieves the minimum. Since the variance 
of f is given by 

A^r = A^p + A^q ICp ^, (5) 

the left-hand-side of Q can also be considered as a function of the three second moments. 

The inequalities Q to (Q and the search for their minima arise from one single math¬ 
ematical problem: 

Does a given smooth function of the second moments have a lower bound? 

If so, which states will saturate the inequality? 

In this paper, we answer these questions for a quantum particle with a single spatial 
degree of freedom by presenting a systematic approach to studying uncertainty relations 
derived from smooth functions f{A^p, A^q, Cpq). Proceeding in three steps we 

1. identify a universal set of states & which can possibly minimize a given functional 
f{A^p,A^q, Cpq)} 

2. spell out conditions which determine the extrema of the functional f as a subset of 
the universal set, (o{f) C (f; if no admissible extrema exist, the functional has no 
lower bound; 

3. determine the set of states ^{f) Q <^{f) which minimize the functional /, leading 
to an uncertainty relation in terms of the second moments. 

The inequalities studied here will be preparational in spirit: they apply to scenarios in 
which the quantum state of the particle \f>) is fixed during the three separate runs of the 
measurements required to determine the numerical values of the second moments. These 
inequalities do not describe the limitations of measuring non-commuting observables 
simultaneously. 
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The paper is divided into two major sections. In Sec.|^we introduce uncertainty func¬ 
tionals and explain how to determine their extrema and minima. The impatient reader 
may jump directly to Sec.|^where we derive new families of uncertainty relations and de¬ 
termine the states minimizing them. We conclude the paper with a summary and discuss 
further applications. 


2 Minimising uncertainty 


To begin, we introduce the uncertainty functional ITTl 


][xp] =/(AVA\Cpg) -A((i/;|t/;) -1), 


( 6 ) 


which sends each element \xp) of the one-particle Hilbert space Jtf to a real number de¬ 
termined by the real differentiable function f{xi, X 2 , xf) of three variables. The Lagrange 
multiplier A ensures the normalization of the states. The variances of position and mo¬ 
mentum are defined by 

= {ip\pfy)) - {^>\p\xp)^, (7) 

etc., and the covariance of position and momentum reads 

Cpq = ^{ip\ ipq + qp) \^p) - {ip\p\ip){tp\q\ip). (8) 

The second moments form the real, symmetric covariance matrix 


_ / Cpcj 

V Cpq 


X w\ 

zv y) ' 


(9) 


with state-dependent matrix elements x = x[-ip) etc. The covariance may take any finite 
real value, re G K, while the variances of position and momentum take (finite) positive 
values only, x,y > 0. States of a quantum particle with vanishing position (or momentum) 
variance and diverging momentum (or position) variance are not taken into account since 
they only arise for non-normalizable states which cannot be prepared experimentally. 
Nevertheless, position (or momentum) eigenstates can be approximated arbitrarily well 
by states within the set we consider. 

It will be convenient to work with states in which the expectation values of both 
momentum and position vanish, {xp\q\tp) = {xp\p\xp) =0. This can be achieved by rigidly 
displacing the observables using the unitary operator 


tfl, = exp [i {poq - qop) /h] , 




( 10 ) 


where po = {xp\p\tp), etc. This transformation leaves invariant the values of the second 
moments (0 and (|^ and has thus no impact on the minimization of the functional / [xp]. 

A lower bound of a functional / [xp] of the form (|^ will result in an uncertainty relation 
associated with the function/(vi,V 2 ,V 3 ). To determine such a bound, we apply a method 
used in lTTHl2j|l3l (see also |[3[l5j[T6l). First, we derive an eigenvalue equation for the 
extrema of the functional J[xp] which also need to satisfy a set of consistency conditions 
given in Sec. 2.3 Then, we introduce a "space of moments" to visualize these results 
(Sec. 


A) and, finally, we determine the minimizing states whenever the functional is 
guaranteed to be bounded from below (see Sec.|2.5|and Sec.|^. 
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2.1 Extrema of uncertainty functionals 

When comparing the values of the functional / at the points \xp) and + e) = \xp) + £\e), 
for any unit vector | e) and a small parameter e, we find to first order that 

}[xp + e]-J[xp]=£D,J[xp]+0{^), (11) 


where the expression 


De = {e 


S 

S{tp\ 


+ 



( 12 ) 


denotes a Gateaux derivative. If the functional }[xp] does not change under this variation. 


= {e\ {x,y,w) - A\xp)'^ +c.c. =0, 

it has an extremum at the state \tp). More explicitly, this condition reads 


/df 3x df 5y df Sw 
\dxS{xp\ dy 3{tp\ dwS{tp\ 



+ C.C. = 0, 


(13) 


(14) 


which should hold for arbitrary variations. Since the vector \e) and its dual (e| can be 
varied independently (just consider their position representations e*{x) and e{x)), the 
expression in round brackets must vanish identically which implies that the complex 
conjugate term will also vanish. Using 


3x _ 3A^p 
S{ip\ ~ S{tp\ 


3{xp\ ^ 


(15) 


a similar relation for Sy/3{tp\, and the identity 


Sw 

S{ip\ 


1 

2 


{qp + pq) \tp), 


(16) 


we arrive at an Euler-Lagrange-type equation. 


The parameter A can be eliminated by multiplying this equation with the bra {xp\ from 
the left and solving for A; substituting the value obtained back into Eq. ( [17) |, one finds a 
nonlinear eigenvector-eigenvalue equation, 

[fxp^ + + y (# + R)) W) = ifxX + fyy + fww) \xp) , (18) 


using the standard shorthand for partial derivatives. 

Eq. is our first result following from the approach conceived in ITT] : the extrema 
of arbitrary smooth functions of the second moments are encoded in an eigenvalue equa¬ 
tion for a Hermitean operator quadratic in position and momentum. However, the equa¬ 
tion is not linear in the state \xp) because the quantities x,y,... ,fw are functions of expec¬ 
tation values taken in the yet unknown state. Previously, similar results had only been 
found for specific uncertainty functionals such as the product of the standard deviations 
or position and momentum. 
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Let us briefly illustrate the crucial features of Eq. ( [l^ in a simple case before systemat¬ 
ically investigating its solutions. For a function linear in v, y, and w, the derivatives fx, fy, 
and fio will be fixed constant numbers. In this case, the operator on the left-hand-side 
of pfll represents a quadratic form in the position and momentum operators, falling into 
one of three possible categories flT l. Up to a multiplicative constant, the operator will be 
unitarily equivalent to the Hamiltonian of (i) a harmonic oscillator with unit mass and 
frequency, (ii) a free particle, p^, or (iii) an inverted harmonic oscillator, p^ — q^. 

In the first case, the spectrum of the operator will be discrete and bounded from below 
(or above); the spectra of the operators in the other two cases are continuous which is 
tantamount to the absence of normalizable eigenstates. Thus, a linear function f{x, y, w) 
possesses a non-trivial bound only if it gives rise to an operator in ( [l8] l which is unitarily 
equivalent to the Hamiltonian of a harmonic oscillator. Generally, our method will signal 
the absence of lower bounds corresponding to the cases (ii) and (iii). 


2.2 Universality 

To find a lower bound of the functional J[ip], we will determine all its extrema and then 
pick those where J[ip] assumes its smallest value. However, Eq. ( |l^ is not a standard 
eigenvalue equation: even for a linear function /, the right-hand-side of ( [l8] l depends 
non-linearly on the as yet unknown state \ip), and if the function / is non-linear, the 
operators on the left-hand-side of the equation acquire state-dependent coefficients given 
by its partial derivatives. 

Nevertheless, the eigenvalue problem can be solved systematically, in a self-consistent 
way. Initially, we treat the expectations v, y, w, and fx, fy, fw in ( |18| | as parameters with 
given values, i.e. we ignore their dependence on the state \tp). The solutions \ip{x,y,w)) 
will depend on these parameters which means that the solutions must be checked for 
consistency since the relations 0 now require that x = {ip{x,y,w)\p^\ip{x,y,w)), etc. It 
may or may not be possible to satisfy these restrictions on the parameters. 

To begin, we write the operator on the left-hand side of ([18| in matrix form. 




fx fw/2 

fw/2 fy 



F z. 


(19) 


Williamson's theorem ITSi ensures that any positive or negative definite matrix can be 
mapped to a diagonal matrix by conjugation with a symplectic matrix E. We will assume 
from now on that the matrix F is positive definite. The negative definite case is easily 
dealt with by considering —f{x, y, w) instead of f{x, y, w). Applied to the 2x2 matrix F, 
Williamson's result states that we can write 


E'^-F 


E = cl, 


E e Sp(2,R), c > 0, 


( 20 ) 


where I is the identity matrix, whenever F > 0 holds. This requires both 

det I = fxfy- fl/4: > 0 and fx> 0, 


( 21 ) 


implying that fy>0 will hold, too. These requirements clearly agree with the observa¬ 
tions made for linear uncertainty functionals f{x, y, w): since the operators p^ and p^ — q^ 


result in matrices F with zero or negative determinant, the left-hand-side of (181 cannot 
be mapped to an oscillator Hamiltonian by means of a symplectic transformation. 

A direct calculation shows that the matrix F is diagonalized by the symplectic matrix 


E = (SyGi,) where 
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G»=(; »), and S,= (7 ° P 2 ) 

with real parameters 

and . = ,23, 

leading to c = \/det F in Eq. ( [^ . The symplectic matrices S-y and G;, give rise to the 
Iwasawa (or ^ JV) decomposition of the matrix G Sp(2,R) (cf. IH^ , for example) 
if they are written in opposite order and the parameter h is replaced by he ^; the third 
factor happens to be the identity. 

Next, we observe that the linear action of the matrices G and S on the canonical pair 
of operators can be implemented by conjugation with suitable unitary operators, 

known as metaplectic operators fl^ . We have, for example. 


1 0 
h 1 

or, in matrix notation. 


P '\ _ ^ibp^/2h f P ^ 
^ J V 4 y 


(24) 


where the unitary operator 


Gfo • z — Gfo z Gt 


Gfc = 

describes a momentum gauge transformation. Similarly, the squeeze operator 


5 ^ _ gh(qp+pq)/2h ^ 


(25) 


(26) 

(27) 


symplectically scales position and momentum according to 


S^-z = S^zS^. 


With L ^ = S^yGfo in (201, we rewrite (19' as 


z"*"- F • z = V detF (S.y ■ G{, • z)"'" • (S.y ■ G{, • z) 


(28) 


(29) 


Finally, using the identities (251 and (28 ' and multiplying Eq. (18 ' with the unitary 
from the left, the condition for the existence of extrema of the functional / [xp] takes on the 
desired form. 


Thus, the solutions 

= S^Gllxp) (31) 

must be proportional to the eigenstates \n),n G No, of a unit oscillator, i.e. a quantum 
mechanical oscillator with unit mass and unit frequency. Equivalently, the candidates for 
states extremizmg the functional }[xp] are given by the family of states, 

\n{b,j)) = GbSy\n), hGNq. (32) 
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Upon rewriting the operator these states are seen to coincide with the squeezed num¬ 
ber states | 6 J. As shown in in Appendix the product of a squeeze transformation S-y 
(with real parameter 7 ) and a momentum gauge transformation Gj, equals 


GhS^^S{^)RU), 


i.e. the product of a rotation in phase space. 


(33) 


= a: e [0,271), (34) 

and a squeeze transformation (with complex ^) along a line with inclination 9, 


S(^) = ^ = re^^ e C. (35) 

Summarizing our findings, we draw two conclusions: 

1. The complete set of solutions of Eq. (j^ coincides with the squeezed number states, 

00 00 

{\n{<x,0) = fj{0\n),<x,^ec} , (36) 

n =0 «=0 


where non-zero expectation values of position and momentum have been reintro¬ 
duced via the translation operator Ta (see Eq. (10 ) and irrelevant constant phases 
have been suppressed. 


2. The value of the right-hand-side of Eq. (301 can take only specific values, 

n e No, 


+ yfy + wfw ^ 1 w ^ 


2 '\/detF 


(37) 


given by the eigenvalues of the unit oscillator. This relation constrains the state- 
dependent quantities of the left-hand-side which needs to be checked for consis¬ 
tency, just as Eq. ( |2T] | does. 

We have thus obtained our second main result. The extrema & of an arbitrary functional 
][tp] characterized by a function f{x,y,w) are necessarily squeezed number states, a set 
which is independent of the function at hand. In other words, the set S’ containing all 
the states which may arise as minima of an uncertainty functional }[ip], is universal. The 
minima of any functional must be a subset S (/) C S which will depend explicitly on the 
function f{x,y,w), determined by the consistency conditions to be studied next. 


2.3 Consistency conditions 


We now spell out the conditions which must be satisfied by the states |n(&, 7 )) in (32 1 - 
or, equivalently, the states |n(a, ^)) in ( j^ - to qualify as extrema for a specific functional 


1 . Recalling that x = A^p, etc., the relations 


and 


V = {n{b,j)\p^\n{b,'y)) , y = {n{b,'Y)\f\n{b,-f)), 


w = -{n{b,j)\ipq + qp) \nib,^)), 


(38) 


(39) 


represent three, generally nonlinear consistency equations between the second mo¬ 
ments since the parameters b and 7 are functions of x, y and w (cf. Eq. (231). 
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2. The values of the moments x,\j and w must satisfy Eq. 

3. The matrix F of the first derivatives must be positive definite. 
Using ( |^ and ( |28| , the first consistency condition in ([38]l leads to 


X = {n{b, j)\p^\n{b, j)) = {n\p^\n) = ^ ^ 

or, recalling the definition of 7 in (231, 

x-\/detF = ?i[/y, neNo- 

Similar calculations result in 


and 


yv^detF ^ (n + - ) ?j/x, neNo, 


— 2 ri 7 '\/detF = fn + ^ j nGNo, 


respectively These conditions may be expressed in matrix form. 


( 40 ) 


(41) 


(42) 

(43) 


F C 
V^detF 



n e No, 


(44) 


involving both the covariance matrix C and F. 

Taking the trace of the last relation shows that Eq. ( [^ is satisfied automatically. With¬ 
out specifying a function f{x,y,w), no conclusions can be drawn about the validity of 
Eqs. (|4T||4^ or the positive definiteness of the matrix F. 


2.4 Geometry of extremal states 

We visualize the interplay of the consistency conditions by expressing them in the form 


= yfy> ^fw = -2wfy, 


and 


XI/ — = ( n + ^ j , 


n e No, 


(45) 


(46) 


following easily from either (41 431 or ( |44| . The third constraint is universal since it does 
not depend on the function/(x,y,a 7 ). Using the variables 

1 1 
w = - (x-hy) > 0 , v = - {x-y) eR, 

we define the three-dimensional space of (second) moments, with coordinates (u,v,w). For 
each non-negative integer, the third condition 


2222 
u — V — w — e„, 


e„= {n + -]fi, n e No, 


(47) 


determines one sheet of a two-sheeted hyperboloid, located in the "upper" half of the 
space of moments, i.e. u > 0 and ? 7 , re G R (cf. Fig.j^. The points on the n-th sheet, which 
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Figure 1: Space of (second) moments, with points (w, v, w): the extremal states of smooth 
functionals J[tp] are located on a discrete set of nested hyperboloids & = USJLo the first 
three of which are shown, using light (n = 0), medium (n = 1) and dark shading (n = 2), 
respectively. The accessible uncertainty region for a quantum particle is given by the points 
on and inside of the convex surface Sq : /4 which coincides with the 

minima of the RS inequality, i.e. squeezed states with minimal uncertainty 


intersects the w-axis at w = +e„, are in one-to-one correspondence with the squeezed 
states originating from the number state |n), forming the set (?„ in 

The consistency conditions ( |4^ clearly depend on the function f{x, y, w) at hand. The 
constraints will only be satisfied for specific subsets (/) of points on the hyperboloids 
Sn, resulting in the /-dependent set of states 

00 

^(/) = U 

n=0 

which contains all the candidates possibly minimizing the functional / [xp]. The candidate 
sets S’{f) may depend on one or two parameters, or contain isolated points only. If the 
consistency conditions cannot be satisfied, then the functional }[tp] has no lower bound. 
Furthermore, if the matrix F is not positive definite for any of the states in the 

method makes no predictions about the minima of the functional / [ip]. We have, however, 
not found any non-trivial cases of this behaviour. 

Finally, we need to evaluate the functional / [xp] for all candidate states & (/) and pick 
the smallest possible value. The states achieving this minimum value constitute the so¬ 
lutions -#(/) C (^(/) of the minimization problem. In their entirety, the minima 
may consist of isolated states or of sets depending on one or two parameters. Usually, 
the states saturating the bound are located on the sheet Sq. 

Let us briefly introduce the concept of a space of moments, defined by the triples of 
numbers (w, v,w)~^ G R^. For n = 0, Eq. ( |47| | is equivalent to ^ which implies that not all 
points of this set can arise as moment triples. The accessible part of the space, bounded 
by the extremal hyperboloid So defined in Eq. ([47||, is called the uncertainty region (cf. 
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Fig.[^. The boundary of an analogously defined uncertainty region for a quantum spin 
s lUoT is not convex. The relation between moment triples and the underlying pure or 
mixed states of quantum particles is discussed elsewhere ETI . 


2.5 Known uncertainty relations 

To illustrate our approach we re-derive three of the four bounds mentioned in the intro¬ 
duction: the uncertainty relations by Robertson-Schrodinger, by Heisenberg-Kennard, and 
the triple-product inequality 

Robertson-Schrodinger uncertainty relation. Defining 

{x, y,w) = xy — up-, (48) 

the matrix of first-order derivatives associated with the quadratic form ([T^ is given by 



and, interestingly, its determinant 

det V = xy — TuP = {x, y, w) (50) 


coincides with the original functional. At this point of the derivation, it is not yet known 
whether the matrix F is strictly positive. 

The relations Eq. ( [4^ do not constraint the parameters x, y, and w, since they are satis¬ 
fied automatically, leaving Eq. ([4^ as the only restriction. Since the left-hand-side of ( [4^ 
coincides with the function f^^Jx,y,w), all squeezed states are candidates to minimize 
the RS functional, 

= (51) 

Therefore, the function comes with the largest possible set of candidates to minimize 
it, given by the union of the sets <§’„ in Fig. The lower bound on now follows directly 
from combining ([4^ with ( [4^ , 



reproducing the RS inequality. Eqs. ( p0| l and ( |52| ) imply that the determinant of F is 
positive ever}rwhere in the uncertainty region. 

The hyperboloid closest to the origin of the {u,v,w)-sTpace provides the states mini¬ 
mizing the function 

(53) 


i.e. the set of squeezed states based on the ground state |0) of a unit oscillator. This is, of 
course, a two-parameter family since the relations take the form 


b = 


z 

/ 

X 


and 7 = - In ( 2 x) 


(54) 


meaning that, with x > 0 and z G IR, both b and 7 take indeed arbitrary real values. 
Thus, each squeezed state can be reached and, when adding phase-space translations, 
we obtain the four-parameter family of all squeezed states as minima of 

= ( 55 ) 
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This property singles out the RS functional among all uncertainty functionals. 

If an uncertainty functional associated to a function f, is different from the RS func¬ 
tional, the first two consistency relations will, in general, not be satisfied automatically 
but impose non-trivial constraints on the second moments. Therefore, the extrema of the 
functional must be a proper subset of those of the RS functional, i.e. S{f) C S, as the 
following example shows. 

Heisenberg's uncertainty relation. Let us determine the minimum of the product of the 
standard deviations Ap and Aq by considering the function 

f^{x,y,w) = (56) 

Its partial derivatives satisfy 




f^ = 0 

J W ^ ' 


resulting in a positive definite diagonal matrix F, namely. 


F 


1 (\/y^ 0 ^ 

2 V 0 ’ 


Eqs. (|^ and ([4^ collapse into the conditions 


detF = 


1 

4 


> 0. 



n G No, 


(57) 


(58) 


(59) 


which determine the value of the product of the standard deviations at the extrema of 
f^(x,y,w), labeled by the positive integers. In the {u,v,w)-space, the intersections of 
the surfaces defined by ( |^ and the hyperboloids ( |47| consist of hyperbolas in the (w, v)- 
plane containing the points (e„,0), n G Nq. The union of these hyperbolas define the set 
S’if^), corresponding to the potential minima of the function f^(x,y) (cf. Fig. §. The 
third condition, Eq. ( [4^ , implies that w = 0. Combining ( [5^ with ( [^ , we obtain the 
bound 

1 \ h 
2)^^ 2' 

reproducing Heisenberg's uncertainty relation Q. 

The family of states minimizing Heisenberg's uncertainty relation is found by using 
the identity (591 for n = 0 in the definition of the parameter 7 in leading to /y = x/h. 


f^(x,y,w) = (^n + 


Since the consistency conditions do not impose any other condition on the variance x, it 
may take any positive value implying that 7 G IR. Since fa; = 0 leads tob = 0, the set of 
states minimizing the Heisenberg's uncertainty relation is given by squeezed states with 
real squeezing parameter. 


^aif) = fccSif^) = {fj^\0), a G C ,7 G ]R} , (61) 

where we have re-introduced arbitrary phase-space displacements. 

Triple product inequality. Using Eq. we see that we need to find the minimum of 
the expression 


f^{x,y,w) = xy (x + y + 2w) 


(62) 
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Figure 2: Hyperbolas in the (w, !7)-plane through the points n G Nq, stemming 

from intersections of the hyperboloids ( [47| and the surfaces defined by ( [5^ . The union 
of the hyperbolas defines the set S (/^) which represents the location of all possible min¬ 
ima of the function the points on the "lowest" (darkest) hyperbola cor¬ 

respond to the set of states ^(/^) which saturate Heisenberg's uncertainty relation. 

in order to reproduce the triple product uncertainty relation (Q. The first consistency 
condition in Eq. ( [4^ implies that x = y; using this identity in the second condition, one 
finds 

x{zv + x){w + ^) = 0. (63) 

Recalling that variances are always non-zero for normalizable states, x > 0, the corre¬ 
lation w must equal either —x or — v/2. According to the first case would imply 
A^r = 0, which is impossible since the operator f has no normalizable eigenstates. There¬ 
fore, using the solution w = — x/2 of ( [^ and x = y in the third consistency condition, 
one finds that ^ 

^ we No, (64) 

must hold. It is now straightforward to evaluate f^{x,y,w) at its extrema to find its 
global minimum. 


f{x,y,w) =x^ = 




which reproduces Q. It is easy to confirm that the matrix F is positive definite with 
determinant detF = / 3. Since the minimum occurs for n = 0, the values of the second 

moments are given by 


fi fi 

X = y = —2w = —= = T- . 
^ ^3 2 


(65) 


These relations fix the values of the parameters in (231, 


b = 


1 

2 


and 



( 66 ) 
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Figure 3: Candidate states = S’^ possibly minimizing the product of three vari¬ 

ances f^{u,v,w), represented by dots located on the intersections of the hyperboloids 
( [47| and the planes defined by the consistency conditions ( |^ ; the point closest to the 
origin, ^(/^) = represents the state |Ho) achieving the minimum of the triple 
product uncertainty relation (|^ (and of any other Ss-invariant inequality associated with 
a functional in (|^). 


Using p2| or (361 we obtain one single state which saturates the triple uncertainty, namely 

|So) = GiS:i,JO)=S,i,3|0). (67) 

If one includes rigid phase-space translations, the set of states minimizing the triple un¬ 
certainty is finally given by the two-parameter family 


,(/^) = {t|Ho),aeC} , 


( 68 ) 


in agreement with IT3]| . 

Geometrically, the state |So) arises from the intersection of the sequence of hyper¬ 
boloids with the surfaces defined by 

1 

M = Te„, v = 0, w = —-Ten, neNo- (69) 

The planes defined by constant values of u have concentric circles in common with the 
hyperboloids, and the vertical ww-plane (given by z? = 0) intersects with each of the 
circles in two points only Finally, the condition on the variable w selects a single one of 
the points with the same value of u. According to ( |^ , the candidate states in {u,v,w)- 
space are located on a straight line. 


^if) 



[ 0 ] ,n e No > , 

V-1/2/ 


and the state |Ho) corresponds to the point closest to the origin (see Fig. §)■ 


( 70 ) 
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3 New uncertainty relations 

3.1 Generalizing known relations 

The linear combination of second moments 


f^{x,y,w) = }ix+ vy+ 2Aw, y,v,A,ElR, 


(71) 


leads to the uncertainty relation 

+ v/t^q + 2XCpq > yv — 


y,v > 0, jiv > 


(72) 


The constraints on the parameters follow from the matrix F in ( [l^ being strictly positive 
definite. The consistency conditions (45' associated with relate both y and w to x 
according to 

y A 


y = -^> 


Then, Eq. simplifies to 


{yv - A2) 


W = -X . 

V 


1 


x^ = I n-h - j , 


(73) 


(74) 


which is consistent due to detF = yv — A^ > 0. Expressing the functional f^[x,y,w) in 
terms of x only, we obtain the bound given in (721, 


f^{x,y,w) = — ---" X = 2enh^J yv — A^ > yv — A^. 


(75) 


Up to phase-space translations Ta/ a single squeezed state saturates the bound, namely 


^(/") = jlF.uA) =G,S^^^^yO)|. (76) 

When expressing the correlation term Cpq in terms of the variance A^r according to 
Eq. (|^, we obtain, for /r = v = 2A = 1 in ( [7^ , the triple sum uncertainty relation 

A^p -6 A^q + A^r > Vsh , (77) 

derived in lil^ , and the minimum is achieved for the state |1,1,1/2) = |Ho) which also 
minimizes the triple product uncertainty (cf. Eq. ^7\ and Eig.|^. 

Sums of powers of position and momentum variances are bounded from below accord¬ 
ing to the inequality 


y{A^pr+v{A^qr > y 


2mm' 
fl \ m-\-m' 


V tfl \ 

}i m 


^ Xi yn in+ni' 

.1/ m' 


?M,m' G N, 
(78) 

reducing to the pair sum uncertainty relation (j^ in the simplest case {y = v = m = m' = 
!)• 


Next, we study a generalized RS-uncertainty functional (481, 

fm%'V'= (w)*" - ’ y>0, m,m' G N. 


(79) 
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For arbitrary integers m and m', the consistency conditions cannot be solved in closed 
form. Setting m' = 2m and assuming that both m > \ and }i > \ hold, we obtain the 
explicit bound 





(80) 


An interesting special case of occurs for m = 1/2 and 0 < pi < 1, 

ApAq-^\CplJ\>^^Jl-pl^, (81) 

which can be treated in spite of the presence of the non-differentiable term. The extremal 
states depend on one free parameter. 


Mfii/i) = \ penn . Jn) \ , V > 0 ; (82) 

in the absence of the covariance term, pi = 0, they reduce to the squeezed number states 
with a real parameter known to extremize Heisenberg's inequality. 

Next, we present an example of an uncertainty relation which seems to be entirely out 
of reach of traditional derivations. Defining the functional 

f{x,y)=x + pie^^^, pi,v>0, (83) 


we obtain the inequality 

> (1 + 2W{h/A^piv)) ^ (84) 

using the fact that Lambert's W-function W(s), defined as the inverse of s( W) = W exp TV, 
is a strictly increasing function. In the limit of p — > 00 and assuming pi = v, the left-hand- 
side of ( [84| ) turns into {pi -T A^p -|- /P^q -T G{\lpi)) while the expansion of its right-hand- 
side produces the correct bound {pi + h + ff{l/pi)), since TV(s) = s -T i^(s^). 

The position and momentum variances at the extremum with label n G No are given 
by 


V = 2pTV 


/ e„h \ w 


enh \ 


and 


e„h \ 

2 ^) ' 


V = 2 vW 

\^V 

respectively. Using Eqs. with detF = {pi/v)e^^'', one finds 

/ enh 


b = 0, 


j = ]\n(^) + Iw 


\2^ 


(85) 


( 86 ) 


(87) 


which means that only a single state (and its rigid displacements) will saturate the in¬ 
equality ( [^ . If p = 1 /, we recover x = y = h^/A as well as & = 7 = 0, i.e. the ground 
state of a unit oscillator since TV(0) =0. 


15 








Let us point out that some general statements can be made about functionals of the 
form / = f{xy,w) and / = /{fix'” + vy’" ,w), i.e. generalizations of the expressions in 
( [78| and ( [7^ , respectively. By examining the consistency conditions one can show that 
the extrema of the first expression come as a one-parameter set, while they are isolated or 
a one-parameter family in the second case. However, without knowing the explicit form 
of the functions no further conclusions can be drawn. 

Rational functions of the variances such us 




^myin 

yx"^' + vy’^'' 


}i,v > 0, 


( 88 ) 


lead to an interesting class of uncertainty relations. If m = m' = 1 / 2 , for example, the 
extremal states of the functional are given by the set 


2,1/2) 


||a,n) = f^S 




(89) 


However, none of these states minimizes the uncertainty functional. In any family of 
squeezed states with real squeezing parameter, the product of the position and momen¬ 
tum variances is constant while their sum increases without a bound. Consequently, the 
functional can be made arbitrarily small without ever reaching the value zero. 


3.2 Uncertainty functionals with permutation symmetries 

The triple product uncertainty relation and the one derived by Heisenberg possess dis¬ 
crete symmetries. Here we investigate more general uncertainty functionals which are 
invariant under the exchange of three and two variances. 


Ss'invariant functionals 

Consider a function of three variables which is invariant under the exchange of any pair, 

/(3) {x, y, z) = (y, V, z) = (v, z, y) . (90) 

We now derive the lower bound of a large class of Ss-invariant uncertainty functionals 
Jlxp] and show that their minima coincide with the state |Ho) minimizing the triple prod¬ 
uct inequality. The variables v, y and z will denote the variances of the operators p, q and 
f, respectively. 

More specifically, we study the minima of sums of completely homogeneous polyno¬ 
mials of degree n, with arbitrary non-negative coefficients, 

N 

p^\x,y,z) = Yj Yj ajke>0, (91) 

n—l j+k+i=n 

dropping the unimportant constant term flooo- The determinant of the associated F- 
matrix. 


det ¥ = fxfy + fyfz + fzfx , (92) 

is positive definite since x, y, z > 0 and the partial derivatives are just positive polynomi¬ 
als. 
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The symmetry under Ss-permutations ( [^ implies that the coefficients must satisfy 
the conditions 


ajke = akji = ajtk, 0 < i,k,e <n, (93) 

so that the first terms of the polynomials are given by 

= flioo(x + i/ + z) 

+fl200 + y^ + + «110 + yz + z:r) 

+«300 + y^ + z^) + fl2io {x^{y + z) + y2(z + v) + z^(x + y)) 

Tfliiixyz + ... (94) 


If the only nonzero coefficients are flioo = 1 or Am = 1, we recover the functionals associ¬ 
ated with the triple sum ( [77| or the triple product inequality Q, respectively In general, 
a completely homogeneous Ss-symmetric polynomial in three variables of degree n > 1 


consists of Kfi = 


{n+3y+6 

12 


terms where the floor function [sj denotes the integer part of 


the number s: each term arises from one way to partition j k + i = n objects into three 
sets with j,k and i elements, respectively [22l. Thus, a symmetric polynomial of degree 


up to N depends on K„j independent coefficients if one ignores the constant term. 

The main result of this section follows from rewriting the consistency conditions 
and (46 1 in terms of the variables x, y and z. 


xfx -yfy + {x-y)fz = o, (95) 

zfz-xfx + {z-x)fy = 0, (96) 

2 {xy + yz + zx) — x^ — y^ — = {2n + l)^h^, (97) 

where we have used the identity z = v -|- y -|- 2zv given in (j^. The conditions ( [^ and 

( [^ imply that the extrema of any symmetric polynomial (x, y, z) occur whenever the 
three variances take the same value. 


X = y = z. (98) 

To show that x = y holds we pick any nonzero term aj]^£x)y^z^ in the expansion ( [^ 
and assume that the powers of x and y are different, i.e. j 7 ^ k; the case j = k will be 
considered later. Due to the symmetry under the exchange v y, the sum also must 
contain the term with = Uj^. Defining t{x,y,z) = Uju [x^y^ + x^y^) z^, the 

first two terms of ( |^ take the form 

xt^ - yty = (;■ - k)ajki [xiy'^ - x’^yi'j z^ . (99) 

Assuming that j = k + S, with ^ > 0, we find 

Xtx - yty = ay+ski^ (x^ - y^) x^y^z^ (100) 

= {x- y)5a^+ski + x^~^y -6 ... -6 xy^^'^ + z^ ( 101 ) 
= {x-y)g+{x,y,z), ( 102 ) 

where g_|_(x,y,z) > 0. Using this expression in ( |^ , the consistency condition takes the 
form 

{x-y){g+{x,y,z) + tz) = 0, (103) 
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with another positive function tz{x,y,z). If < 0 we write k = j — 5 = /+|^| and 
eliminate k instead of j from (991, only to find that its left-hand-side again turns into 
{x — y) multiplied with a positive function. If the powers of x and y of the term a^^x^y^z^ 
are equal, j = k, one immediately finds that {xd^ — ydy)apxiyh^ = 0, also reducing 
Eq. (|^ to (x — y)dzajkix>yh^ = 0. 

The argument just given covers all terms in the sum ( |M] |, and the positivity of the 
coefficients implies that the first consistency condition can only be satisfied for x = y. 


(S') 

Using the symmetry of f^'{x,y,z) under the exchange y ^ z, an identical argument 
leads to the identity y = z. 

Using ([^ to evaluate the left-hand-side of Eq. ^7) results in 


X = y = z = Tenh, n G No, (104) 

where t = \/473, so that we obtain an uncertainty relation for any Ss-invariant function 


f^^\x,y,z) > f^^\x,y,z) 


x=y=z=ti / \/^ 


(105) 


This result correctly reproduces the special cases of Eqs. (Q and ( [77| , and there is only 
one state which saturates the inequality, namely |Eo) given in Eq. ^7) . Letting N —> oo 
in Eq. ( |M] |, we conclude that the main result of this section, Eq. (|105|), also applies to any 

(S') 

Ss-symmetric function (x, y,z) with a Taylor expansion with positive coefficients and 
infinite radius of convergence, as long as its first partial derivatives exist. 


S2-invariant functionals 

Assume now that, in analogy to Eq. ( [^ , we have a functional depending on just two 
variances in a symmetric way. 


/®(r,y)=/ra(y,^). 


(106) 


An argument similar to the one given for the function (x, y, z) results in the uncer¬ 
tainty relation 




x=y=h/2 


(107) 


which covers the cases of Heisenberg's relation ([TJ and the pair sum inequality (|^. Thus, 

) 


(2\ 

the actual form of the function at hand determines whether the set of minima 
will depend on a continuous parameter or not. If the functional is invariant under scal¬ 
ing transformation x —)■ Ax, y —)■ y/A, in addition to the permutation symmetry, there 
is a one-parameter family of solutions and the right-hand-side of Eq. ( 107| | achieves its 
minimum on the set of points with xy = {h/2)^, not just those with x = y = h/2. 

To derive (107||, we suppose that the function f^\x,y) has an expansion in analogy 


to f^\x,y,z) in Eq. ([^ but without the variable z. Adapting the reasoning applied to 


fjy \x,y,z), the consistency equations ( |45| are found to imply x = y and w = 0. Using 
this result in ( |4^ , the bound x > follows immediately, so that the inequality ( |107| | 
must hold for S 2 -invariant functionals. 
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4 Summary and discussion 

This paper responds to the fact that, without exception, the states minimizing the known 
preparational uncertainty relations for a quantum particle in one dimension are given by 
(sets of) squeezed states. To explain this fact we systematically study lower bounds for 
smooth functions /(A^p, A^q, Cpq) of second moments. The resulting theory explains the 
universal role of squeezed states for preparational uncertainty relations depending on 
second moments only, and it completely charts the landscape of inequalities of this type. 

The chain of inclusions 

3 ^{f) 3 M(/) (108) 

concisely summarizes the general structure of our findings. First, we have shown that 
only squeezed number states of a quantum mechanical harmonic oscillator with unit 
frequency and mass occur as extrema of an uncertainty functional }[xp] depending on 
second moments. We denote this universal set of states by S’. Second, the extrema of a 
specific functional Jf[ip], associated with a function/(A^p, A'^q,Cpq), form the subset S{f) 
of the universal set S. Third, the functional will assume its minimum for one or more of 
the extrema S‘{f), a subset which we denote by ./#(/). The set of minima may be empty, 
.y#(/) = 0. If it is not empty, a lower bound on the functional }f[ip] has been found, and 
it represents a preparational uncertainty relation in terms of the second moments. 

Strictly speaking, we obtained the relations Sc^ ^ '^a(/) 3 instead of Eq. ( |108| l 

as it is possible to move quantum states in phase space without affecting the values of the 
second moments. The four-parameter set Su = TaS, for example, consists of the squeezed 
states S plus those obtained from them by means of the translation operator defined 
in Eq. ®. Thus, each state saturating a specific inequality with vanishing expectation 
values gives rise to a two-parameter family of minima. 

Our results have a useful geometric representation in the real three-dimensional space 
of moments. The uncertainty region, consisting of all triples of moments which can arise 
from (pure or mixed) states of a quantum particle, turns out to be a convex set bounded 
by a one-sheeted hyperboloid. The boundary is invariant under the elliptic rotations, 
hyperbolic boosts and parabolic transformations which generate the group SO(l, 2). This 
observation squares with the importance of the group SO (1,2) in quantum optics where 
coherent and squeezed states are ubiquitous. 

The invariance of the bounding hyperboloid can, in turn, also be understood as an 
invariance of the functional }[ip] defining the surface. Each point on this hyperboloid is 
associated with a unique Gaussian state saturating the Robertson-Schrodinger inequality. 
Changing from an active view of transformations (i.e. mapping one state with minimal 
uncertainty to another one) to a passive view, we see that the RS-functional (or a suitable 
smooth function thereof) is invariant under the elements of the group Sp (2, R) applied 
to the canonical pair {p,q). The allowed symplectic transformations include rotations, 
scalings and linear gauge transformations, or shears. 

Repeatedly, we have encountered subsets of the maximal possible symmetry group 
Sp(2,R) (cf. IZl). The Heisenberg functional f^{x,y,z) in Eq. ( |5^ , for example, is in¬ 
variant under a scaling transformation p ^ Xp, q ^ q/X, with A > 0, resulting in a one- 
parameter set of states with minimal uncertainty depicted in Fig. Similarly, uncertainty 
functionals invariant under permutations of order two or three are minimized by states 
with corresponding symmetry properties. Examples are the triple product uncertainty 
relation Q and, more generally, the functionals with discrete symmetries discussed in 
Sec. IQ 

We also derived new and explicit uncertainty relations. Fig.j^uses the (&, 7 )-plane to 
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Figure 4: States on the boundary of the uncertainty region minimizing known and new 
uncertainty relations parameterized by the real numbers {b, 7), with h = 1; each point of 
the plane corresponds to a squeezed state saturating the RS-inequality (|^; points on the 
vertical dashed line represent minima of Heisenberg's uncertainty relation Q; the two 
curved dashed lines indicate the minima of the modified RS-inequality ( [sT] ) with m = 1/2 
and values }i = 1/2 (bottom) and = 1/10 (top); the full dots correspond to minima of 
S 2 -invariant functionals ( |107| | such as the pair sum (|^ and Ss-invariant functionals ( 105| l 
such as the triple product Q. 


illustrate the sets of states which minimize (some of) the uncertainty relations discussed 
in this paper. The sets of minima ^ may depend on two parameters (all squeezed states 
minimizing the RS-inequality), on one parameter (such as the real squeezed states satu¬ 
rating Heisenberg's uncertainty relation) or consist of a single point only (associated with 
Ss-invariant inequalities such as the triple product inequality, for example). We have not 
been able to backward-engineer functionals which would be minimized by prescribed 
subsets of the plane such as a circle or a disk. The minimizing states we found were all 
pure, located on the boundary of the uncertainty region. In principle, functionals could 
also take their minima inside this region although we have found only trivial examples 
with this property l 211 . 

We currently investigate three conceptually interesting generalizations of our ap¬ 
proach. First, there is no fundamental reason to restrict oneself to uncertainty functionals 
depending only on second moments in position and momentum Il23l . On the contrary, 
higher order expectation values would enable us to move away from Gaussian quantum 
mechanics which is largely reproducible in terms of a classical model "with an epistemic 
restriction" of the allowed probability distributions Il24l . Including fourth-order terms 
{tp\q^\tp), for example, will result in an eigenvalue equation ( [18| which is not related to a 
unit oscillator in a simple way. It is known that fourth-order moments for single-particle 
expectations can give rise to inequalities which cannot be reproduced by models based on 
classical probabilities Il25ll26ll . Thus, it might become possible to study truly non-classical 
behaviour in a systematic manner using suitable uncertainty functionals. 

Secondly, our approach can be generalized to the case of two or more continuous 
variables. We expect that a systematic study of uncertainty functionals becomes possible, 
leading to criteria which would detect pure entangled states. For example, a generaliza¬ 
tion of the triple uncertainty relation Q to bipartite systems has been shown to detect 
entangled states in a quantum optical setting 127]. Other scenarios are known which rely 
on intuitive choices of suitable bi-linear observables ||2^|29| . 
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Finally, the comprehensive study 1201 of uncertainty relations for a single spin s has 
been limited to observables transforming covariantly under the group SU(2). The method 
proposed here is easily adapted to investigate functionals depending on arbitrary func¬ 
tions of moments. 
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A A Baker-Campbell-Hausdorff identity 


The relation Gi,Sy = S{^)R{x) in Eq. ( [^ can be shown by requiring that both products 
map the annihilation operator a = {q ip) /to the same operator. We obtain 


A A.J. q. / _ U 

GbSryaSyGfj = a ( coshy — i-e 


^sinhy + i^e'^ 


(109) 


and 


S{^)R{x)aR^{x)S'^{^) = ae ^^coshr — d^e '-'^e'^sirihr, ( 110 ) 

respectively. Equating the coefficients of the operators d and d^ leads to two equations 


h 

cosh 7 — = e cosh r , 

sinhy -|- sinhr. 


( 111 ) 

( 112 ) 


which we need to solve for the variables ^ and x- Separating the real and imaginary 
parts of the first equation, one finds that 


X = arctan 


i- ^ 

l + e-^'T j ^ 


(113) 


In a similar way, the second equation allows one to solve for the function tan(0 — x) 
which, upon using ( 113| l, leads to 


6 = arctan 


( ^ \ { ^ \ ^ 

- -^ -I- arctan --^ G (— tt , ti . 


1 - e-27 

The case of 7 = 0 needs to be treated separately leading to the relation 


71 


0 = ±—-|-arctan e( — tt, tt). 

Finally, the condition cosh r cos x — cosh 7 results in the expression 

02 ^ 1/2 


(114) 


(115) 


r = arcosh 


/ 

( cosh^ 7 -|- j S [ 0 ^ 00 ) / 


(116) 
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which establishes the desired identity 

The number states \n),n G N, are eigenstates of phase-space rotations R{x)- There¬ 
fore, the product G{,S-y acts on those states according to 

GhS^\n) ^ S{0\n), (117) 

where an irrelevant phase has been suppressed. Thus, the operator S(^) generates all 
squeezed states from |0) when the parameter ^ runs through the points of the complex 
plane. 
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